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1. Introduction 


Formulating string/M theory in a time-dependent background remains an elusive prob¬ 
lem. The only observables in string theory are S-matrix elements, this is certainly true 
in the perturbative formulation of string theory where conformal symmetry on the world- 
sheet plays the role of the guiding principle in constructing consistent asymptotically flat 
background, it is also true in a nonperturbative formulation of M theory, the matrix the¬ 
ory, where scattering amplitudes among DO-branes and their bound states are assumed to 
exist. However, S-matrix does not exist for most of interesting cosmological backgrounds, 
it certainly does not exist for our universe. Perhaps, a reformulation of observables is the 
key to extending string/M theory to include time-dependent backgrounds. 

The matrix model proposed by Craps et al. is an attempt to formulate string theory 
in a time-dependent background [|l| , the metric in this model depends on a null coordinate 
and in the Einstein frame it exhibits a null singularity at the “big bang” point. This model 
was subsequently generalized to a class of more general backgrounds in [^, and to a class 
of even more general backgrounds in [^, |^], and 0 (a concrete model in this class was 
previously studied in detail in @). For related work on time-dependent backgrounds, see 
0, i, i, and 0. 


So far, except for the decoupling argument presented in 0 , there has been no indepen¬ 
dent check on the correctness of the matrix proposal. The effective action of a DO-brane 


in the background generated by another DO-brane was derived in [|n|, where it is noticed 
that the usual double expansion in the relative velocity v and the inverse of the relative 
separation b fails when time is sufficiently close to the big bang point. Although there is no 
dehnition of scattering amplitude between two DO-branes too, we believed that it makes 
sense to talk about the effective action at later times. In the present work we shall make 
the usual one loop calculation to see whether we can obtain the small velocity expansion 
of 101 . To our surprise, we shall see that the w^ term in the one-loop calculation does not 
vanish and is complex. This is a rather astonishing result. 

We are faced with two possibilities, our result may indicate that the matrix proposal 
is incorrect, or it may signal an instability of the two DO-brane system at later times, since 
the f ^ term in the effective action is complex. However, we can not locate a physical reason 
for this instability at present. 


The layout of this paper is as follows. We use the background held method of [|0] to 
write down a gauge-hxed action and expand it to the second order. We compute the one- 
loop contribution of the off-diagonal huctuations to the effective action of two DO-branes 
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in sect.3 when the relative velocity vanishes, and hnd it eqnal to zero. The term in the 
one-loop contribntion is calcnlated in sect.4, and we hnd a non-vanishing complex term. 
We show that the small velocity expansion makes sense in the hat matrix theory and the 
terms cancel in appendix A. Appendixes B is devoted to discnssions on the propagators. 


2. Basic setup 

In 1^], the anthors consider a hat type IIA backgronnd, with a nnll linear dilaton, 
(j) = —Qx^. The Einstein metric has a cnrvature singnlarity at = —oo. A matrix string 
action is proposed in [|l|] to describe the theory nonpertnrbatively. The type IIA backgronnd 
can be obtained by compactifying M theory on a circle, along the ninth direction. In 
the backgronnd is lifted to M theory, and the corresponding matrix theory is BESS like 
T^ . DO brane interaction is fonnd out by considering the shock wave solution in fH 


To get the shock wave solution, the authors have compactihed the ninth direction and 
averaged the source over that direction. The Routhian of a graviton in the presence of 


another is 


\p- p'_v'^r where p_ is the null momentum of the 


test particle, {2%)“^RR' p'_ is that of the source particle, and Cn is some hxed numerical 
coefficient, especially, ci = 1, C2 = ■^. R is the radius of the M-theory and R' is that of 
9th direction. From the form of the Routhian, one can see that there is no static potential 
between two gravitons, and there is no correction, either. In this approach, we are 
expanding the effective action in terms of p'_v‘^r~^. Note that 




2 2 


271-i?' 


is just the physical gravitational constant in IIA string theory. Therefore it is clear 
that the expansion is a supergravity perturbation. 

In the present paper, we will just consider the case of two DO branes, and hence 
P- — and p'_ = ^ 2 -nR 2 -KR' • shall in this paper work in the scheme of ||], and 
use the matrix model action instead of the matrix string action. We compute the effective 
potential of two DO branes with separation both in the ninth direction and in the transverse 


directions. To compare our matrix model calculation with the supergravity result in |n 
the separation in the ninth direction should be integrated out in the end. 

The matrix theory action includes the bosonic part Sb and fermionic part Sp- Set 
the Planck scale /p to 1, the two parts can be written as 




W]^}, 

4 2 


Sf= dtTriie'^DtO - Re^^e'^'y^[X\ 9] - Re^'ygiX^, 6]}, 


( 2 . 1 ) 
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where i,j = 1, • • - 8, runs over the eight transverse directions, and Dt = dt -\- i[A, is the 
covariant derivative. Rescale t Q —> and to absorb 

the R in the action, we have 


Sb = j dtTr{{DtX^f + e-^^\DtX^f + ^e^^^[X\X^f + }, 

Sf= [dtTr{i9'^Dt9 - e^^9^'y,[X\9] - 9^-fg[X^,9]}. 


( 2 . 2 ) 


To calculate the effective potential, we use the background field method [13|. Expand 


the action ( p.2|) around the classical background field by setting X^ = + Y^, 

jj. — 1,2 ■■ - 9. The fiuctuation part of the action is a sum of five terms 


S 


S, + S9 + Sa + Sf ermi T Sghost- 


(2.3) 


In the following, we will determine the explicit form of each term, 
choose the gauge 


G = dtA- ie^^^[B\ X^] - i[B^, = 0. 


It is convenient to 
(2.4) 


In the standard gauge fixing procedure, we need to insert 

(2.5) 

into the path integral, where ^ is a gauge parameter, f(t) is chosen to be f{t) = for 
later convenience, g{t) is any function. The path integral is independent of the choice of 
g{t), so we can multiply the path integral by J[dg{t)]e~'^^^^^ . A/p is given by the variation 
of G under gauge transformation, independent of g{t). Thus by changing the order of 
integration, we can integrate out g{t), and get a gauge fixing term 

Sgf = -e-^^^G^. ( 2 . 6 ) 

So the bosonic action of the fiuctuation is 

Syi = J dtTr{{dtY^f + e‘^^\[B\ Y^f + Y^f + [B\ Yf 

+ 2[B\Y^][Y\Y^] + ^[Y\Y^f)}, 

Sy9= j dte-^^^^TriidtY^f + e^^\[B\Y^f + e-^^^[B^,Y^f 

+ 2[B\Y^][Y\Y^]+2[Y\B\Y\Y^] + [Y\Y^f)}, (2.7) 

Sa= J dtTr{-e-‘^^\dtAf - e-^^^[A,B^f - [A,Bf 

+ 4idtB^[A, y*] + Aie~^^^dtB^[A, - AQie-^^^B^[A, 

+ 2idtY^[A,Y^] - 2[A, B^][A,Y^] - [A, Y^f 
+ 2ie-^^^dtY^[A,Y^] -2e-‘^^^[A, B^][A, -e“2Q*[A, T^]^}. 
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Since we are considering two DO-branes, the Yang-Mills fields are jnst 2x2 matrix. 
We will choose the backgronnd to be diagonal, 






( 2 . 8 ) 


The backgronnd for A and other transverse directions are chosen to be zero. This cor¬ 
responds to, in comoving coordinate, two zero-branes moving towards each other with 
relative velocity v in the direction, and transverse separation b in the direction and 
c in the x^ direction. Write the matrix in terms of U{2) generators. 


W = i(Yo*l2 + WV“), + YaV“), 

A = ^e^^(Aol2 + A^a^, 6 = ^(0ol2 + ^aCT“). 


(2.9) 


where a = 1, 2, 3. The 0 components in this decomposition describe the free motion of the 
center of mass and will not be written explicitly in the following. Then np to qnadratic 
terms, the actions for the finctnations are 


Si = ^J dt{Yi\-dt^ - - c^)Yi^ 

+Y2\-dt‘^ - b^e^^^ - - c^)Y2^ 

+Y3^(-dt^)Y3^}, 

S9 = ^j dt{Y^\-dt^ + _ b2^2Qt _ ^2^2^2Qt _ ^2)y-^9 

+Y2\-dt^ + - c^)Y2^ 

+Y3\-dt^+Q^)Y3^}, 

^A = ~ J dt{Ai{-dt'^ + _ c‘^)Ai 

+A2{-dt'^ + <5^ - _ c‘^)A2 

+A3i-dt^+Q^)A3 

+4ve^\AiY2^ - A2 Yi^) - 4Qc{AiY2^ - 


Define 


r? 




( 2 . 10 ) 


( 2 . 11 ) 
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Then the actions for A and Xg become 
Sa+ = ^ Jdt{A'^i{-dt‘^ + Q^ - - i2Qc)A'^i 

+A'^2{-dt‘^ + Q^ - - i2Qc)A'^2 

+A3i-dt^ + Q^)A3 + i2V2ve^^[{A+i - ^“ 1 )^ 2 ^ + (^+2 - 
Sa- = ^ j dt{A-i{-dt^ + Q^- + i2Qc)A-i 

+A~2{-dt'^ + + t2(5c)A“2 

+Y3^{-dt^ +Q^)Y3^}. 

( 2 . 12 ) 

Dehne new fermionic helds, 

^+ = ^(^i+t^ 2 ), e- = ^{0i-i02). (2.13) 

Then the action is 

Sf = J dt6'^{idt + nte^* 7 i + 66*5*72 + 679 ) 6 '+ + ^ 03 ^{idt) 03 . (2.14) 

The ghost action is determined by the inhnitesimal gange transformation of G, 

Sg = j dtCi*{-d^ - - c^)Ci + C2*{-d^ - b^e^^^ - - c^)C2 

+ C3*{-df)C3. 

(2.15) 

Before doing any calcnlation, we can see that the flnctnation action for is inde¬ 
pendent of the separation, and hence has nothing to do with the interaction of the two 

zero branes. We will not consider them in the following. 

3. Static case 

First we will analyze the sitnation when n = 0. This corresponds to two zero-branes 
static in the comoving coordinates. To calcnlate the one loop interaction, we need to 
integrate ont the qnadratic flnctnation, which can be written in the form of determinants, 
det“^(-a2-6262<5*-c2), for y+ 2 ^ t = 1 , ■ ■ ■ 8 , 

det“5(_a2 g 2 _ ^ 2 g 2 Qt _ ^2 _ ^ 2 Qc), for A+i,A+ 2 , 

det “2 (—52 _|_ q 2 _ ^2g2Qt _ ^2 ^2Qc), for A~i,A~ 2 , (3.1) 

det(—— 6^6^*5i _ ^ 2 ^^ I'qj, Ci^ 2 , 

det {idt + 66 * 5*72 + 679 ), for 9+. 
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We use Schwinger proper time formalism to calculate the determinants. For any Hermitian 
operator A, the determinant is represented by 

roo 


/■“ds 

d = ln(det A) = - / —Tre- 
Jo s 


(3.2) 


(3.3) 


Thus we need to calculate the heat kernel, K{t'^t\s) =< t'\e >. K{t,t'\s) satishes 

the differential equation and the boundary condition, 

idsK{t', t] s) = AK{t', f; s), 

K{t', t; 0) = 5{t — t'). 

For the hrst determinant in (|3.1|) , A = —d^ — — (?. To solve (|3.3|) , we hrst solve the 

static shrodinger equation 

= {-dl - - c^)y\{t). (3.4) 

The two linearly dependent solutions of (|3.4|) are Bessel functions 

J±k(x), for K ^ Z, or Jk(x), Y^{x), for k e Z, (3.5) 

where x = = A + c^. Yn{x) has singlarity at x = 0, and are not in 

consideration. Since the operator A is hermitian, A is real, and k is either real or pure 
imaginary. 

Using an integral of Bessel function (eq. 6.574.2 of jlSfl ) 

(3.6) 


'dx T , ^ 2sin7r(V^) 

X Tl[IJ, + — U) 


an orthonormal basis can be constructed, 
Vo^it) = 


Qu) 


['7iti;(x) T T—j(^(x)], cj ^ 0, 


2 sinh(7rc<; 

fn{t) = \/ 4 :QnJ 2 n{x), n = l,2, 

To check the orthogonality, 

QoJ 


(3.7) 




u ^sinh[|(w + w')] 


2 sinh( 7 rct;)yo Qx 

e 


\Jiuj-\-e{x) T T—(x)] [T—(x) T (x)] 

sinh[|(w - w')] e 


+ 


sinh(7rw)' {oj + u’) 7r[(^)2 + e2] {u - u') 7r[(^)2 + e2] 

= 5{uj — uj') + d(ci; + uj') 

= 5{uj — uj'), 

nOO 

/ dtfn{t)fm{t) = Sn,m, 


} 


' —oo 
poo 


dtff^{t)fn{t) (X sin(?T.7r) = 0. 


(3.8) 


6 













We have deformed izw by a small real part, iiw —iiw + e and used the identity d{z) — 
limg^o —nV^T- In the fourth line w > 0 is taken into account. 

To check the completeness, we will need to prove that all Jk{x), > 0, k ^ 2n,n G 
Z_|_ can be expanded in the basis. Dehne 


J^{u) = / dtyl{t)J^{x), J” = / dtfn{t)J^{x), 

Jo Jo 

_ «O0 C)0 

J^{x) = / dujJ^{uj)y^{t) + ^ Jnfn{t)- 


n=l 


(3.9) 


Using (|3.6|) , one hnds that 

-T-4cosh(^) sin(^) 




Hence, 


2(5sinh(7rci;) 7r{K'^ + uj'^) 


duJ^{u)y^{t) = / du- 


= 4, 


Ur(-l)"sin(7rf) 
Q ti{k^ — A'n?) 


sin(^)wJi^+e(a:) 


JO J-oo 7rsinh(^)(K2+ cj2)' 

The large order behavior of the Bessel function is 


(3.10) 


(3.11) 


J^(x) 


(3.12) 


Then the integral (P.ll[) can be evaluated by closing the contour in the lower half plane. 
Simple poles are aX u — —2m, —in. 


duJ^{u)y^{t) = J^{x) -sin( —) ^^(-1)’ 

TT Z 

n=l 
oo 

= W - 


n 




— Ati? 


J2n{x) 


(3.13) 


n=l 


Therefore, 

Ji^{x) = Jk(x). (3.14) 

When K = Aiiuj + e, the above equations still hold. Then the other linear combination of 


J±uj{,x), Jiu}{x) — J-i^{x) can be also expanded in terms of the basis (|3.7|) , and so are not 
included in the basis. In fact, we have shown that any normalizable eigenfunction can be 
expanded in terms of this basis, which is enough to guarantee that the basis is complete. 
(The completeness of this set of the eigenfunctions was discussed previously in [pT| .) 
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Then the heat kernel can be expanded in terms the orthonormal basis, 


K{t\t;s) = / duy^{t')*yu;{t)e ^ 

3.1. The bosonic effective potential 


,i[{2Qnf+c‘^]s 


(3.15) 


Having fonnd ont the heat kernel, it is straight forward to write down the determinant 
explicitly, by 5i = — ^ Jff^dtK{t,t\ s). The trace in (|3.2|) is now an integral over t. 

In order to compare with the result obtained on the supergravity side [|^, we need to 
compactify the 9-direction and smear the result over the circle. This is equivalent to sum 
the images in the covering space and then average over the compactihed circle. On the 
matrix theory side, we need to calculate the one loop effective potential of two DO-branes 
separated also by c in the direction, integrate over c and then divide by 2tiR'. This 
procedure is expected to give us the result that is to be compared with our earlier result in 
0 . Here R' is the radius of X^. Now there are altogether four integrals in our calculation 


of determinant, the integral over t, s, uj, and c. We can first do the the integral over c. Then 
the smeared determinant becomes 


5.= 


-1 


2\/—'KiR'jQ 


r 

2\/—'KiR' Jo 


ds 

si , 


ds 

si , 


dt{ 


duj 


Qoj 


’ —oo 
oo 


2 sinh(7r(3; 




+ '^4QnJ^^{x)t 


dt{ 


Qoj 


n=l 
poo 

duj 

2smh(7r(3;) 

OO 


i{2Qr 




x[l - i{Qijo) s H-]4- 7 ^4QnJ2^(x)[l i{2Qn) s H-]}. 

n=l 

(3.16) 

Here we have extended the integral range of u from (0, cxo) to (—cx),cx)). We use the 
notation = e^ 2 “, and ln(z) = We have rewritten the exponential in 

the form of power series. Using the large order behavior of Bessel function, we have 


U) 


2n+l 


sinh(7r(3; 




U> 


271-|-1 


sinh(7r(3;) 


exp[— 2 z( 3 ;lnz — In ( 3 ;] = 




to 


2n+l 


sinh(7r(3;) 


JLix) 


^2n-\-l ^ 

——-r exp[2fc<; -|- 2zc<; In — 

smh(7r(3;) 2 

exp[2i(3; -|- 2iu) In | — 2iu) In)^; -|- ttu] 

sinh(7r(3;) 


sinh(7r(3;) ’ 

2iuj In ( 3 ; — 2iui In i — In uA 


(3.17) 
















Close the contour in the lower half plane, we can see that the integral of each term pro¬ 
portional to at the inhnity is zero. Because of the third line of (|3.17|) , we will 

meet a divergence at inhnity in each term proportional to Jiuj{x)J-iuj{x). Note that this 
divergence is independent of the x, and therefore can be subtracted. Then 


oo Quj 

du) ~ : ; Z Jiuj{^X^\^Ji^{x^ “1“ T—S “(“ * * *], 


2 sinh(7rct;) 

= -4QnJ2n{x)[l + i{2Qn)^s + 


(3.18) 


The above just cancels with the summation in ( b.l6|) term by term. Although we are not 
sure about the convergence of the expansion, the exact cancelation of each term between 
the integral (|3.18|) and the summation in (|3.16|) has show that 5^ = 0 up to a physical 
irrelevant constant. So the bosons coming from the i directions give no contribution to the 
effective potential. 


For the second and the third determinants in (|3.1j). The heat kernel becomes 


K+{t\ t-s)^ / ^ fn{t')fn{t)e 


i[{2Qn)^ 


n=l 
oo 

2i„ X—^ , V I / ■r^\2^ 


K_{t', t-s)= I b + ^ +(c-.Q) 

Jo 

(3.19) 

Then take the same procedure as in eqs. ( b.l6| ), and (|3.18|) , we will hnd that the 

bosons coming from the gauge held and give no contribution to the effective potential, 
either. 

The ghost determinant is the same with that of A*, and hence give the same result 
except for a minus sign. 

In a word, we hnd that there is no static potential coming from the bosons. 


3.2. The fermionic effective potential 

In the fermionic sector, there are 16 degrees of freedom for each SU{2) index. Since 
there are only three gamma matrix relevant here, we can choose a basis to make the gamma 
matrix and the held block diagonal, 

7i = 0-2 ® Is, 72 = CTs ® Is, 79 = 0^1 Is- (3.20) 


Dehne 

K^p{t\t-,s) =< t'\exp{-iAf{i)s)\t >a(3, 
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where Aj = idt + be^^'f 2 + C 79 , a, f3 = 1,2 label the two 8 x 8 block matrix. Then 
Kajsit'it, s) satishes the following differential equation and initial condition 


idsKaf3{t',t]s) = {Af{i))c,pKpf3{t',t;s), ,t]0) = - t). 


(3.21) 


To hnd the solution, we write 

K 2 (}{t',t\ s) 


dXKi 3 {t', t] A)e 


—iXs 


In the following, we just write Kp for short. Then from (|3.21| ), Kp satisfy the following 
differential equations 

(idt + - X)Kip + cK2p - 0, 

(3 22) 

{idt - be^^ - X)K 2 P + cKip = 0. 


These equations are equivalent to 

cK2p^-{tdt + be^^-X)K^p, 

{d^ + b^e^^^ + + 2iXdt - iQbe^*)Kip = 0 . 


(3.23) 


Denote X/{Q) by u, and c/Q by c'. Take the ansatz Kip = f{w^c\t')'Ll:{oj^c\t). Then 
/(w, c', t') factorize and the equation for '(/^(w, c', t) has two linearly independent solutions, 
a;(“^/^-*‘^)Mi/ 2 ,ic'(~ 2 fa:) and Where Mx,p,{x) is a Whittaker 

function. A general solution for (|3.23|) is 




+ 




(3.24) 


where f{LJ,c',t') and g{uj,c',t') are chosen to satisfy the initial condition. Then 


J-00 ^ —Atyxx 


— \M 1/9 ■ >(-2ix') ( ^V 2 ,ic'( 2za;) A 




J-00 ^ —Aix/xx' 


A[Ml, 2 ,-.A-‘ 2 ^x') I 


Mii2,ic'{-‘^iA A 
M-i/2,ic'{-‘^ix)) 
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(3.25) 









When s = 0, w can be integrated ont and gives d{t — t'). Then 


where 


0 ) 


Sc,p5{t-t') 


D{x, x', c') 
—A.i\/xx' 


D{x,x',c') = Mi/ 2 ,ic'{-‘^ix)M_i/ 2 -ic' 


{-2ix')+Mi/2-ic'{-2ix)M_i/2^ic'{-‘2ix'). (3.26) 


In appendix B, we will prove that D{x, x, c') 
in ( CT) . 

Finally, 


—4ix. So (|3.25| ) satishes the initial condition 


df = - 


dc 


ds 


— OO 

poo 


27tR' n S 


dt 


'0 

dc 


trKa,,p{t,t'; s) 

roods roo 


'—OO J—oo‘^^^ Jo 


djjjQ exp[—zQws] 


^ —OO 

X [Mi/ 2,ic'(-2ia;)M_i/2_ic'(-2ix) + ^i/2,-ic'(-2*a;)M_i/2,ic'(-2*a;)] 

n (it n 


(3.27) 


' —OO J —OO - 


'0 


dwQ exp[—zQws] 


The tr in the hrst line means a trace of the 2x2 matrix and an integral over t. The integral 
is divergent but is independent of b and a. In fact, this is just the phase shift generated 
by a free operator Regularize the phase shift by subtracting — Hq is 

the Lagrange for free fermions. We can see that the fermions give no contribution to the 
effective potential. 

Then we can draw the conclusion that there is no static effective potential. 


4. Effective interaction at the order 

Now we are going to investigate the case when there is a small relative velocity be¬ 
tween the zero branes. Since it is difficult to compute the determinants directly, we will 
perturbatively expand around n = 0. 

From (|2.10|) , ( pj.l2|) , (|2.14| )and (|2.15|) , we can see that the only possible terms term 
odd in V in the perturbation series come from the the fermionic action (|2.14|) . These terms 
vanish because the trace of odd number of Gamma matrix is zero.We shall in this section 
calculate various terms. 


11 


















Denote the terms coming from the hrst order bosonic contribntion by 6 i, the second 
order bosonic contribntion by 62 , the hrst order ghost contribntion by ( 71 , and the second 
order fermionic contribntion by f 2 - Then, 




+ {At{t)A+{t)) + {A+{t)A+{t)) + {A-{t)A-{t)) + {A-{t)A-{t))], 

/ OO 1 poo 

—J dtoW«*i(c,(t)cr(f)> + {C2(f)Q(i))i, 

1 rir r°° 

^2 =77 / 7r~QJ / / dt2{vV2fe^^^e^^^ 

^J — oo^^^ J— 00 J— 00 

X {[< 1)24"*“ 1(^2) > + < A 2(^1)^ 2(^2) >] < T 2 ^(ti)l 2 ^(^ 2 ) > 

+ [< ^ 1 (^ 2 ) > + < ^"'~2(tl)^"'~2(t2) >] < bi^(ti)Yi^(t2) >}, 

°° dc 


h-\ 


^CXD poo 

L^2nR'J.A'J.^ 


dt2vHit2e^^^e^*^ 


X Tr[7iGy (ti, t2)7iG^/(^2, ^i)]- 


(4.1) 


Since we are only interested in the effective potential, we do not have to do all the 
integrals in the second order contribntions. Dehne ti = t + ^r, t 2 — t — ^r, integrate ont 
r and c, we are left with an integral of t, which combined with the hrst order pertnrbation, 
will give the ehective potential to v'^ order. In Appendix A, we will show how this procednre 
is carried ont when the backgronnd is hat. We hope that this procednre also goes throngh 
here, as we shall see, there is a problem arising at this order. 

The propagator < Y^{t 2 )Y^(ti) >= Gi(t 2 ,ti), satishes the diherential eqnation 


(-7, - - c^)Gi(f2, *i) = iHh - h), 


(4.2) 


and is related to the heat kernel by Gi{t 2 Ai) = ~jQ^dsK{t 2 Ai'i^)- 


Gi{t2,ti) = - dsK{t2,ti,s 
Jo 

. r, Quj 

=l / dOJ -—-r 

2smh(7rw, 

00 


"Ao; (371) [t/—(372) T (3)2)] / [(Qc*^) c te] 


+ i ^ 4 ( 5 nJ 2 n( 3 :i) J 2 n( 3 : 2 )/[-( 2 ( 5 n)^ - - ie] 


n=l 


= - 9{ti - t2) 


TT 


2Q sinh( 7 rc 


A—ic'(^2 ) [ Aic'(^1) T A—ic'(^l)] 


TT 

9(t2 ti) —r : -pr J—ic' ) [A^c' (^ 72 ) T J—ic' ( 372 )]- 

siriiiiTrc ) 


(4.3) 
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From the second line of (14.31) to the forth line, we have integrated w by contour integral, and 
assumed c' > 0. When d < 0, just replace d with —d. Using the asymptotical behavior of 
Bessel function at large order ( p.l2| ), we see that -iujixd) ~ 

When — ^2 > 0, we should close the contour in the upper half plane, and otherwise the 
lower half plane. The poles at ±2ni cancels with the sum in the third line of ( [4.3|) . So 
only the poles at ±(c + fe) contribute to the propagator. It is difficult to obtain a compact 
result. We will take a 6 —oo limit to obtain the asymptotic behavior. Or equivalently, we 
let hj (Q) to be of order 1, and let t > cxd . 




cos(x 2 + \i'nd — cos(xi — j) 
Qdx\X2 sinh(^) 


+ {ti t2). 


(4.4) 


We use Giitd) to calculate 6i, which becomes 


Giitd) ~ 


[sin(2a;) + 1] coth(^) 

2 Qx 


i cos(2x) 

2 Qx 


Dehne 

G+(l,2) =< A+i{ti)A+i{t2) >=< A+2{ti)A+2{t2) >, 
G-{ 1 , 2 ) —< A i(fi)A 1(^2) >=< ^ 2 {ti)A 2(0) > ■ 
They satisfy the following differential equations 

(- 9 j - - (c + iQf)G+{t2, ti) = i 5 {ti - t2), 

- [c-iQf)G-{t2di) = i5{ti -t2). 


The solution G+ can also be obtained from the heat kernel. 


G+{t2di) = - dsK+{t2di,s) 


— iuj {xi)][Ji^{x2) + J—iuj (^ 2 )] 


r 


2 sinh(7rci;) [(Qw)^ ~ (c + iQ)" 
4QnJ2n(Xl)J2n(x2) 


-(2Qn)2-(c + zQ)2 

^(ti ^ 2 )^ 7 - /\rT / \ T 

I T i(c'+i) (^l)J 


2Qsinh[7r(c' + z)] 

“ 2 Q sinh[7r(c'V z)] (^ 2 ) + (a; 2 )] 

cos[x 2 + |z7r(c' + z)-f]cos(a:i-f) 

- O /TWUcosh(^) -+ 


(4.5) 


(4.6) 


(4.7) 


(4.8) 
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and 


G'+(t, t) 


[sin( 2 a:) + 1 ] tanh(^) icos( 2 a;) 

2Qx 


2Qx 


To get G -, just replace c' by —c' in . 
Then the hrst order contribution adds up to 


(4.9) 


-iV, = 


dc 


,271R 


t) + G+{t, t) + G-{t, t)] 


+ sin(2a:)] f°° Qdc' cos(2a:) Qdc' 

-^- i -^- i IHF 

f°°Qdc' , /Ttc'. 
b In ttR' ^ 2 ^ 


(4.10) 


In the last step, we have omitted the trigonometric functions because they are periodic 
functions and fluctuate violently at large argument. The integral of c' seems to give a 
divergence, but this is just caused by our using the asymptotic expansion of the bessel 
function. That step hides the depression of the large cb In fact, from the large order 
behavior ( 3.17|) , we can see that there is indeed no divergence in the cb Hereafter, we can 
just put the this divergence aside. 

To compute 62 , we again need to take the limit |a:i| S> 1 and 1 x 2 1 3> 1. This is 
equivalent to ti 3> 1 and t 2 S> 1. In this case, we use the asymptotic expansion of Bessel 
function when xi and X 2 are large in (|4.3|) and ( fl. 8 |) . Multiply ( |4.3|) and (|4.8|) , and integrate 
out r = ti — t 2 , we will get the dependence on x = t = \iti + t 2 )■ The relevant 

integral is 

-iVb2 = 


dc 


) ttR' 

»oo 


2x^e^^* / dTGi{tiR2)[G+{t2Ri) + G-{t2Ri)] 


ttR' 


dr 


0 


0{ti - t2)v‘^ 
62 


, TTC 


{i coth(-—) [cos( 2 x 2 ) + - sin( 2 xi + 2 x 2 ) 


H— sin( 2 xi — 2 x 2 )] — [1 — sin( 2 x 2 ) + sin( 2 xi) 

2 

- ^ cos( 2 xi - 2 x 2 ) + ^ cos( 2 xi + 2 x 2 )]} + {ti ^ 62 ) 


(4.11) 


' dc' .. , /TTcb T /. \i 

^p-0^coth( —)[/o(4x) - Lo(4x)] - 


Qdr + 2ibo(4x)}. 


Considering xi 3> 1, and X 2 3> 1, The trigonometric functions with argument Xi, X 2 , and 
2 xi + 2 x 2 fluctuate quickly in the x S> 1 limit, they average to zero and hence can be 
omitted. The term J^^Qdr seems to be divergent. However, remember the limit we are 
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taking here, ti, t 2 3> 1, and t fixed, so the range of both ti and t 2 is proportional to t. Then 
the range of r = —12 is also proportional to t, and the term proportional to j^^Qdr is 

finite and increases with t. 

To calcnlate / 2 , we will need the fermionic propagator, defined by GQ:,/ 3 (ti, ^ 2 ) = 
{T 6 j^{ti) 6 '^{t 2 ))oLp- It satisfies the following differential eqnation, 

[idt^ + 66 * 5*172 + c79]G«,^(ti, 62 ) = -iSa,i3S{ti - 62 )- (4.12) 

The propagator is related to the heat kernel Ka/ 3 {t 2 ,ti-, s) ronghly by Ga,/ 3 (ti,t 2 ) = 
dsKai 3 {t 2 :ti] s). Bnt there is some snbtly in determining the time ordering in each 
term. This is related to the bonndary conditions. We are not going to solve the problem 
in this way. Instead, we take the 6—^0 limit. The limiting case will be the propagator for 
massive fermions, which will be analyzed in Appendix A. When 6—^0, the argnment of 
the Whittaker fnnction is small. We will have 

(4.13) 


We determine the time ordered propagator by comparing its small 6 limit with the propa¬ 
gator of a massive fermion in Appendix A , conseqnently 


^ 11 ( 61 , 62 ) = {4:iy/xiX2)~^[9{ti - t 2 )M_i/ 2 ,-ic'{-‘ 2 ix 2 )Mi/ 2 ,ic'{- 2 ixi) 

- 6 (t 2 - ti)M_ii 2 ,ic'{-‘ 2 ix 2 )Mii 2 -ic'{-‘iixi)] 

<^ 22 (^ 1 , ^ 2 ) = {Aiy/xiX2)~^[9{ti -t2)Mii2-ic'{-‘2‘ix2)M_ii2,ic'{-‘iixi) 

-0{t2 - ti)Mi/2,ic'(-2ia:2)M_i/2 _ic'(-2ta:i)] 

(4.14) 

Gl2{tl,t2) = {Aiy/XiX2) ^[9{ti -t2)Mi/2-ic'{-‘2iX2)Mi/2,ic'{-‘2iXi) 

+0{t2 — ti)Mi/2,ic'(~2ta:2)Mi/2,-ic'(~2ta:i)] 

G2liti,t2) = {4iy/XiX2)~^[0{ti -t2)M_i/2-ic'{-‘2iX2)M_i/2,ic'{-‘2iXi) 

+0{t2 — ti)M_i/2^ic'{ — ‘2‘iX2)M_i/2,-ic'{~‘^'^Xi)]. 

Again let b/{Q) to be of order 1, and take the ti ,2 3> 1 limit also, we will get a finite 
integral with respect to both r. In the following, we will need to nse the asymptotical 
expansion of Whittaker fnnction at large argnment. 


M 

M 


i ±ic'(-2*a:) ~ exp(T7rc* + ix)^/- 2 ix 


-i,±ic'(-2*a:) 


r(±tc' -I-1) 
r(±2tc' +1) 
' r(±tc' +1) 


ex 


p(—tx) V—2fir. 


(4,15) 
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Then 


Tr[7iG'/(ti, t2)7iG^/(^25 ^i)] — [G'ii(t25 ^i)G22(^i, ^ 2 ) + ^ 22 (^ 2 , ti)Gii{ti, t2) 

— Gi2{tl,t2)Gi2{t2,ti) — 6*21 (tl, t2)G2l{t2, ^l)] 

= {^XiX 2 )~^e{ti - t 2 )[M‘i^/^ _-^,{- 2 iX 2 )Mli^ -^,{- 2 iXi) 

+ ^i%,-ic'(“2ix2)M^;|^/2,ic'(“2ixi)] + (ti ^ t2) 

cosh(2zx2 — 2ixi + 27rc') 

rs_/ — - 

cosh^ (ttc') 

(4.16) 


The fermionic contribntion to the effective potential is thns 


-iVf^ 


dc 


b T^R'J-oo 


dT^v‘^e^^\t^ - —)Tr[7iG/(fi,t2)7iG'/(^2,ti)] 


d,dc' t) 2 g 2 gt cosh( 27 rc') 


'kR 


li 


TT 


c»hVc') 

i2TTv‘^t^e^^^ sinh(27rc') 


(4.17) 


cosh^(7rc') 


■[/o(4x) - Lo(4x)]} 


From (|4.10|) , ( |4.11|) , and (|4.17|) we see that the effective potential proportional to 
does not vanish. The late time potential contains both a real part and an imaginary part. 
In the following, when we nse “proportional to”, we mean that we ignore some nnmerical 
coefficient, inclnding the the integral of c'. The leading real part comes from the bosons, 
( [4.1U|) , proportional to increases with t. The leading imaginary part also 

comes from the bosonic part, ( [1.11|) , proportional to -ip- J^^Qdr. This term is hnite 
and increases as t as we have explained following eq. (|4.11| ). We may also pay attention 
to the snbleading terms, which are hnite, and may have some physical signihcance. The 


snbleading imaginary part comes from the fermionic contribntion, which is proportional 

9 9 904 / \ 

to —ivte^Ko{Ax) ~ ™ ^ ^ . The snbleading real contribntion comes also 


v 


from the fermionic contribntion, which is proportional to v‘^t‘^e^^^[Io{4x) — Lo(4x)] 




,Qt 




Both the real part and the imaginary part of the effective potential are 


V 

proportional to positive power of Q, so when Q 
potential also vanish as t —> cxo. 


0, both vanish. The snbleading effective 
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5. Conclusion and Discussions 


We study the effective potential between two DO-branes in a time-dependent matrix 
theory at the one loop level. When the two DO-branes have no relative motion in the 
comoving coordinates, we find that there is no effective potential. This result is expected 
if there is supersymmetry, thanks to the cancelation between bosons and fermions. What is 
surprising is that there is no supersymmetry in our case. The bosonic and fermionic phase 
shifts are both divergent but do not dependent on the physical parameter, the seperation 
h. So upon suitable regularization, they are both zero. 

When we consider the case when u 7 ^ 0, the exact form of the effective potential is 
not calculated because of the integrand is too complicated. Too see that the potential is 
non-trivial, we examined the behavior of the potential in later times. The corrections 
do not cancel in one loop calculation. Moreover, there exists an imaginary part in addition 
to a real part. This result seems to contradict with our supergravity calculation. When we 
compactify the X'^ direction, we get a type IIA string theory with string coupling constant 
Qs = and the effective 10 dimensional gravitational constant is oc 

Supergravity loop expansion is in terms of gravitational constant. But we see no sign 
of this expansion in matrix calculation. Furthermore, the imaginary part of the effective 
potential may imply an instability of the 2 DO-brane system. As the two DO-branes move 
apart in the comoving coordinates, certain modes in the two DO-brane system become 
tachyonic, and the imaginary part just signals creation of these modes. 
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Appendix A. Perturbation in flat background 


When the background is flat, BFSS matrix model has been tested to two loops. Here 
we will use our perturbation method to repeat the result to one loop order. Set 6 = 0 and 
Q = 0, we just return to the situation investigated by |^. The u = 0 case is similar. The 
determinants we are going to compute becomes 


det^°(—— c^) 

for 


■■9 and A12, 


det“^(—— c^) 

for 

Cl,2, 


(A.l) 

det“^(fdf - 1 - C79) 

for 

e+. 
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The propagators are Gh(t,t') — —* I for all the bosons and the ghosts. For the 
fermions, 


/\ I 


Gii(t, t') — G 22 {t-> t') — —^{t — t')e 
^12(^0 =^21(^0 = 

where e(t — t') — \ [0{t — t') — 9{t' — t)]. For the v ^ 0 case, we need 

h = -idt-vh\Yi^{t)Yr{t)) + (iT(t)iT(t)) 

+ ("4i(t)Ai(t)) + {A 2 {t)A 2 {t))], 


(A.2) 


gi^il dtvh^[{Giit)GUt)) + {G2it)GUt))i 


' —oo 

/*oo /»oo 


b2 = - 


dtidt2v'^{{Ai{ti)Ai{t2)){Y2^{ti)Y2^{t2)) 


— oo^ —oo 


+ {A 2 ih)A 2 it 2 )){Y^\G)Y^\t 2 ))}, 

r*OG poo 


f2 = 


1 


dtidt2V^tit2 X 8Tr[7iG/(ti, t2)7iG/(^2, ^i)]- 


— oqo —oo 


(A.3) 


In order to get the effective action, we do not need to perform all the integrals. Dehne 
t = T = ti — t 2 , integrate ont r, and snm over all terms above, we will get the 

effective potential before the smearing: 



(A.4) 


The varions factors comes from the connting of degree of freedom. They snm np to zero. 
So there is no term in the effective action. 


Appendix B. The proof of an identity 

Here we will give the proof of the following identity 


D{c, z) = Mi/ 2 ,ic{^)M_i^ 2 ,-ic{z) + Mii 2 ^-ic{z)M_ii 2 ^ic{z) — 2z, 


(B.l) 
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where z is pure imaginary. In the following, we will treat D{c, z) as a function of c, and 
view z as a parameter. Using the steepest descendent method, we can get the large |c| 
behavior of Whittaker function, Mi ^^^^{z) ~ ~ when |c| — 00 . 

So hm|c|^oo-D(c, 2 ) ~ 2z. Using M\^^[z) — — A + |,2// + 1 ; 2 ), and the 

relation $(a, 7 ; z) — e^$(a — 7 , 7 ; —2), we can write iA(c, z) in terms of $ as 


D{c, z) — z[^{ic, 2ic + 1; z)^{—ic, —2ic + 1; —z) + $(ic, 2ic + 1; —z)^{—ic, —2ic + 1; z)]. 

(B.2) 

$(a, 7 ; z) as a function of 7 has single poles at 7 = —n, and analytic elsewhere. Near the 
pole, limjic+i —„*(ic, 2 ic+ l;e) ~ „i(7+i+n) f n + 2 ; «), where 

( ^ - 1) - 1] • ■ ■ |i(« - 1) - «)/[(« + I)!]- 

If n is odd, the above is zero, so the potential poles in the upper half plane are at 2ic+ 1 = 
—2n. However, 


lim D(c,z) 

2ic+l^-2n 

^ / i(2n - 1) 

(2n)!(2m+1+ 2n) V 2n + 1 

x[^2n+l ^_x)2n+l| ^ q_ 


■|.(?7U. 2n + 2; 2" + 2; -z) 


(B.3) 


The same phenomenon happens when 2ic + 1 — 2n. Thus D{c, z) is analytic in 
the complex plane as a function of c. Since D{c,z) approaches to 2z as the |c| —> cx), 
D{c, z) = 2z by Cauthy integral formula. 
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